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Abstract 

Lattice formulation of a fermionic field theory defined on a ran- 
domly triangulated compact manifold is discussed, with em- 
phasis on the topological problem of defining spin structures 
on the manifold. An explicit construction is presented for the 
two-dimensional case and its relation with the Ising model is 
discussed. Furthermore, an exact realization of the Kramers- 
Wannier duality for the two-dimensional Ising model on the mani- 
fold is considered. The global properties of the field are discussed. 
The importance of the GSO projection is stressed. This projec- 
tion has to be performed for the duality to hold. 

Introduction 

The massless Majorana free fermion theory belongs to the same universality 
class as the critical Ising model on a regular lattice [|1], 0, |3|, |J. An explicit 
construction of the Majorana-Dirac- Wilson fermion field theory on a ran- 
domly triangulated plane was introduced in |J. This theory was shown to 
be equivalent to the Ising model also outside the critical region. In ref. || 
Cartesian coordinates were assigned to the nodes of the lattice. The direc- 
tions of the links and of the related gamma matrices were expressed in the 
global frame of the plane. This approach works for lattices embedded in a 
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flat background where one has at one's disposal a global frame of the under- 
lying geometry @, 0, §]• However, if one wants to generalize it to a lattice 
on a curved background where no global frame exists, a field of local frames 



10 , |TT|1 has to be introduced. This being done, one can put fermions on a 



curved manifold with any topology and one can eventually attack, for exam- 
ple, problems of field theory on a dynamical geometry like those encountered 



in string theory or in quantum gravity |L_4| |15|, [16|, |17|, |1<; 

This generalization was partially carried out in 110, |TTJ where an ex- 
plicit construction of the Majorana-Dirac- Wilson operators on curved com- 
pact two-dimensional lattices was introduced. 

Here we extend these studies. In particular, we discuss the significance 
of the GSO projection, which as in string theory also here plays an impor- 



tant physical role fll2 , 13[] . We show that with a careful treatment of the 



global properties of the Dirac operator and of the spin structures on the 
manifold one can find a strict mathematical one-to-one equivalence between 
the partition function of the Majorana- Wilson fermions and that of the Ising 
model. We show explicitly that in our discretization of the Dirac operator 
on a compact manifold, the GSO projection - the summation over all spin 
structures - does remove the non-contractible fermionic loops, that is those 
not corresponding to the domain-walls of the corresponding Ising model. 

Further, we show that for the duality to hold exactly as a one-to-one map 
between the Ising model on a triangulation and on its dual lattice, a sort of 
GSO projection has also to be done. Different spin structures for the Ising 
field are simulated by physical cuts produced by the introduction of antifer- 
romagnetic loops, which mimic antiperiodic fermionic boundary conditions. 

The paper is organized as follows. In section |1] we give an introduction 
to the problem of defining the Dirac operator on a compact manifold. It is 



text-book material [[13|, pO . We recall it here for completeness, to keep the 
article self-contained. In section [2|, we show how to adapt the standard Wil- 
son discretization scheme of fermions on the regular translationally invariant 
hypercubic lattice to the local-frame description, which can be general- 
ized to the case of irregular curved lattices. In section [| using as an example 
the standard toroidal regular lattice, we discuss the sign problem and the 
global properties of the fermionic field on a compact manifold. In section 
[| we argue that in the case of irregular lattices the local frame description 
is particularly natural, and then in section |5| we show how to lift this con- 
struction to the spinorial representation. In doing this we introduce rotation 
matrices between neighboring frames which are crucial for the construction. 
In particular, using the spinorial representation of these matrices we are able 
to define in section ^| the Dirac- Wilson operator. The standard definition of 
the partition function representing quantum amplitudes is recalled in section 
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[7|. In this section we also list the properties of the mathematical expressions 
encountered in calculating the partition function. In section |] we calculate 
the partition function using the hopping parameter expansion. The topolog- 
ical loop sign problem emerges naturally there. The issue of loop signs is 
discussed in more detail in section where the sign is defined as a function 
of classes of loop homotopies. The relation between signs of non- contract ible 
fermionic loops and of domain-walls in Ising model and the topological as- 
pect of the duality is discussed in section jTUJ. In section [□] we give two 



analytic examples, calculating the critical temperature of the Ising model on 
the honeycomb lattice and the critical value of the hopping parameter on 
the dynamical triangulation, making use of the existence of the exact map 
between the Ising model and the fermionic model. We close with a short 
discussion. 



1 Preliminaries 

The aim of this paper is to discretize a theory of fermions on a random, 
possibly fluctuating geometry. Let us first recall some basic facts about the 
continuum formulation of this problem. 

Consider a D-dimensional compact Riemannian manifold, on which a 
coordinate system £ M is defined. If a nonsingular change of coordinates £ M — > 
£ M is performed at some point x on the manifold, then a linear transformation 
of the components of any vector or tensor field in the tangent space at x 
has also to be carried out, in order to ensure the invariance of the theory 
under coordinate transformations. For vectors, the matrix of this linear 
transformation reads : 




Since the change of coordinates is not singular, the determinant of A is 
nonzero. The matrices A thus form a linear group of non-singular real matri- 
ces GL(D, R). The basic difficulty in any attempt to apply the transforma- 
tion law (|T|) to a fermionic field is that the group GL(D, R) has no spinorial 
representation. In other words, one cannot directly apply the information 
encoded in A to transform a spinor when changing the coordinates. In order 
to overcome this difficulty one has to restrict somehow the group GL(D, R) 
to its SO(D) subgroup, which does have spinorial half- integer representa- 
tions. One can do this by introducing an additional field of local orthonormal 
frames. More precisely, at each point x of the manifold one introduces a ba- 
sis e a (x), a = 1, . . . , D, in the tangent space, which obeys e a (x) ■ eb(x) = 5 a b 
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(orthonormality) and ei(x) A e2(x) ... A ec(x) > (orientability), where the 
symbols • and A denote the internal and external products. 

Expressed in a given coordinate system the orthonormality and ori- 
entability conditions read : 

g^(x) e»(x) e u b (x) = 5 ab , e(x) = det e£(x) = yfgjxj > . (2) 

The matrix e%(x) is called the vielbein. It is non-singular, and one can denote 
its inverse matrix by e°(x). Thus one has, for instance, e a ^x)e b v (x)5 ab = 

With these vectors one can also associate gamma matrices 7 a , {7 a , 7 b } = 
25 ab , that can be chosen so as to have the same numerical values 7 a for all 
points x. One can write the Dirac matrices in the curved coordinates as 

7 M (x) = e^(x)7 a . 

The price to pay for introducing this new field is that one also has to 
introduce an additional connection on top of the Levy-Civita connection. The 
new connection u (which is called the spin connection) allows one to calculate 
covariant derivatives of objects that have frame indices. For instance, the 
covariant derivative of the vielbein itself is given by 

VX = ^ + re A e a A -^X (3) 

The reward is that the spin connection can be lifted to the spinorial repre- 
sentation, and we can calculate the covariant derivatives of spinors as well : 

= + ^ ab a ab 4> , (4) 

where a ab = ^[7 a , 7 b ] is the rotation generator in the spinorial representation. 
The action for fermions coupled to gravity can now be written as : 

S= ±/d^e^7"V^ = 1 1 Jd D x^( 1 a -V a )^ 

(5) 

= ±jd D xd D y ip{x) D(x,y)ip{y) . 

The Dirac operator on the manifold is 

D(x,y) = 6(x-y)r(x)-V a (x), (6) 

or, less formally, just 7 • D. We shall discretize this operator in the next 
section. Before doing so, however, let us discuss its topological properties. 

Locally, one can always define a continuously varying field of frames. 
However, doing this globally for a compact manifold is usually impossible. 
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What can be done instead in this case is to cover the manifold with open 
patches, in each of which one can separately define a continuous field of 
frames, and for any region of overlapping patches U and V provide transition 
matrices for recalculating the frames when going from one patch to the other : 

\eu\aix) = [Ruvt [evh(x) . (7) 

Here, the transition function Ruv is a SO(D) rotation matrix. It follows 
that the spinors in the overlapping region can be recalculated as : 

foM«0*0 = [Kw\i foM*(aO • (8) 

where TZuv is an image of Ruv m the spinorial representation. In a region 
where three patches U, V, W intersect, the transition matrices must obviously 
fulfill the following self-consistency equations : 

RuvRvwRwu — 1 , TtuvTtvwTLwu = 1 • (9) 

The second equation can be almost automatically deduced from the first 
one by rewriting it in the spinorial representation. However, because the 
spinorial representation R — > ±1Z is two-valued, the signs of the 7£'s are not 
automatically fixed by i?'s. In other words, one has to adjust in addition 
the signs of the transition functions for the spinors in such a way that the 
consistency equation is fulfilled in any triple intersecting patch. 

This is a global topological problem. If it is solvable on the entire mani- 
fold, the manifold is said to admit a spin structure. In two and tree dimen- 
sions, the question of the existence of a spin structure reduces simply to the 
manifold orientability; in higher dimensions the problem is more complex. 

Another important question is: how many non-equivalent spin structures 
are admitted on a given manifold ? In two dimensions, the answer is 
where g is the genus of the manifold This number is related to the 

number of possible sign choices for independent non-contract ible loops on 
the manifold. 

A good discretization scheme should reflect all these topological proper- 
ties. As will be seen, the explicit construction for two-dimensional compact 
manifolds to be proposed in the present paper does fulfill this requirement. 

The Dirac operator (||) can be expressed in local coordinates as 7 A *V A1 , or 
alternatively in frame components as 7 a V a , i. e. without reference to local 
coordinates. The construction proposed in this paper is, in fact, coordinate- 
free : we shall express everything in frame indices a, without referring to 
coordinate indices \i. 

In the lattice construction, the nearest neighbor relation that mimics the 
structure of the continuum formulation will be given by a local vector : at 
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each point i on the dual lattice we shall define local vectors nji pointing to 
the three neighboring vertices j. To calculate derivatives (differences) in the 
direction of we shall decompose it in the local frame ei a . Similarly, all 
vector, tensor and spinor indices of objects from the tangent spaces will be 
expressed in these local orthonormal frames. Lifting the construction from 
the vector to the spinor representation of the rotation group, we shall store 
the information about nearest neighbors in the form of rotation matrices. We 
refer to them as to the 'basic rotations', and denote them by the letter B. 
The advantage of using rotations is that we can express them in the spinorial 
representation, B — > B. 



2 The discretization scheme 

Let us start with a discussion of fermions on a regular flat lattice, using 



the Wilson formulation ||22|| . Then, we shall see how to go over, after some 



modifications, to the case of irregular lattices. 

The Dirac- Wilson action for free fermions reads 



S = -y E { W 1 + + *s(l - V 4 )***} + \ E ( 10 ) 



t./i 



where the multi-index i describes the node position on the lattice, and Jl is 
one of the D directions of the lattice. The gamma matrices 7 M are rigidly 
associated with these directions : 

{7*, 7"} = 2*"". (11) 

In the Euclidean sector, the Dirac field is represented by independent Grass- 
mann variables ^> a and a — 1, . . . , N. In particular, for D = 2, the 
dimension of the spinor representation is iV = 2. In the following, spinor 
indices will usually be implicit; we shall write them explicitly only when 
necessary. 

We shall now rewrite the action ([K]) in a coordinate-free form which can 
be extended to the case of irregular lattices. 

Instead of using the multi-index i to describe the vertex position, we 
associate with each vertex a single label, say i, which is a coordinate-free 
concept. Obviously, the particular choice of a label does not have any physical 
meaning and the theory has to be invariant under relabelings. The physical 
information will be encoded in the nearest neighbor relations. 

Using these labels, the action can be cast into the following form : 

S = -KY, ViH^j + ~ E > ( 12 ) 
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Figure 1: A hypercubic lattice with translational symmetry and a global 
frame that fixes the coordinate directions for the entire lattice. Alterna- 
tively, one can use local frames that vary from point to point. This has the 
advantage of being generalizable to a curved background. 



where the first sum runs over oriented links connecting nearest neighbors on 
the lattice. The hopping operator Hij is defined as 

H l3 = l -(l+n irl ) } (13) 

where local vector pointing from j to i, being assumed that the two 

are nearest neighbors. Note that in the sum over oriented links, each link (ij) 
appears twice, once as (ij) and once as (ji); since we clearly have = —riji, 
we see that the action (|T2|) is indeed equivalent to (|10|). 

Even at this stage it is more elegant to stop referring to coordinates and 
instead use components of the global frame E a = (E\, E 2 ). Thus, we replace 
7 M by 7 a , and decompose the nearest neighbor vector into components in 
this frame. The product n^- • 7 can then be expressed as : 

Uij ■ 7 = n iiiQ 7 a = n^-17 1 + n ij)2 l 2 . (14) 

Written in the form (0), the action is now coordinate-free, but it still 
depends on the global frame through the vector components and the 
gamma matrices 7°. Such a global frame and a common spinorial basis exist 
only in exceptional geometries, like the regular torus or plane. In order to 
define a theory on another topology or, generally, on a curved background, 
we have to get rid of this concept and use local frames instead. 
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One can introduce independent orthonormal frames as in fig. [I]. At each 
lattice point i one has a pair of orthonormal vectors (e^e^)- In particular, 
on a torus the local frames ej a can be obtained from the global frame E a by 
local rotations : 

e ia =[R t ] b a E b . (15) 

The spinor components are transformed by these rotations into their com- 
ponents in the local bases ipi : 

where the matrices TZi belong to the half-integer representation of the rota- 
tions Ri : 

In component-free notation the equations ([TBI), (HI) and dl7|) read : 

ei = RiE, <4) t = K t m t) $ i = %Tli 1 , HiiRi 1 = R i 'y. (18) 

Using this notation, one should remember that the matrix TZ acts on the 
spinor indices whereas R acts on the frame indices. Using the local frames, 
we can write the action (TO) as : 



S = -K^iUi^ + ^£V#* (19) 

<ij> < 

where 

Hi, = TZiH^Uj 1 = -lit [1 + riij ■ 7] KJ 1 TZ(lZf . (20) 

Here, is a matrix allowing to recalculate the components of a spinor going 
from a frame j to the frame i. In other words, it is a sort of a connection 
matrix that performs a parallel transport of spinors between neighboring 
vertices. 

So far, equation ( p0|) is written in a hybrid notation, because the spinors 
are already expressed in the local frames whereas and 7 are still written 
in the global frame E. However, applying (|17]) to ( p0|) one finds : 

7^ ny • 7 = Kcplli 1 = n ijja R a b 7 fe = ng ] ■ 7 (21) 

where in the local basis the vector has the components 

41 = n ^aRt , (22) 



different from the global frame components . The new bracketed index 
(i) now differentiates between different local frames where the components 
of the vector are calculated; thus, nfj refers to the same vector as , 
but with components expressed in a different frame. Intuitively, what the 
equation means is simply that the components of a vector in a rotated basis 
can be alternatively calculated by performing the inverse rotation on the 
vector itself while keeping the basis fixed. 

An important point is that the crossover from the global description to 



the local one as in (21) preserves the numerical values of the 7 a matrices. In 
other words, 7 1 associated with the local direction en at a point i has the 
same numerical value as 7 1 associated with the e^i at any other point j, and 
likewise for r y 2 . 

Using the components nfJ of the nearest neighbor vector in the local 
frame i, we can now write (|20|) as 

' 'l + rfl-y] U l3 . (23) 



Hij 2 



Alternatively, using the features of nfj discussed above, we can cast the 
hopping operator into several equivalent forms : 



Hij 2 



l + njy-7 = ~ l-nJV-7 U ij = U^- l + ng'- 7 • (24) 



1 



2 



1 



2 



These different expressions for Hij correspond to different ways of calculating 
the hopping term ipfHijipj in fllpp. One method is to first parallel transport 
the spinor ipj from j to i, getting Uijtpj, and then to calculate the correspond- 
ing scalar in the frame i, as is done on the left hand side of (p4|) . Sometimes 
it is convenient to replace = — riji in order to change the direction of the 
vector between indices % and j, as is done in the second expression. Alterna- 
tively, one can first transport the spinor ^ from % to j , which gives ^>M%ji 
and then calculate the corresponding scalar in the frame j, as is done on the 
right hand side, etc. All these expressions are equivalent and can be deduced 
from each other, so that the most convenient one is always chosen. 

The additional upper index in the brackets makes formulae visually less 
transparent but removes the logical ambiguity which otherwise might lead 
to confusion. We will therefore extend this notation to all objects occurring 
in our construction. For example, ipj = Uijipj means that the spinor ipj is 

transported from j to %. Similarly, = ipf^Uij means that ipi is transported 
from % to j. There is no summation over the repeated indices. The only 
exception will be made for objects calculated in the frame belonging to the 
point where they are themselves defined, since in this case leaving out the 
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upper index does not cause any ambiguity. For example, we will write ipi 
instead of tp, 



(0 



Using this notation, the Wilson action becomes : 
S = 



(25) 



Contrary to (10), this form of the Wilson action can now be generalized to 
any random irregular lattice. It also makes direct contact with the continuum 
formalism (|5|). Finally, note that it is invariant under a change of the local 
frames : 



( , -> lUei , % -> , tf« -> , £4, _> KMijKf ■ (26) 



where i?j are arbitrary local rotations, and 7^ are the corresponding matrices 
in the spinorial representation. 



3 A topological problem 

Let us return to the consequences of the fact that the (spinorial) half-integer 
representation of the rotation group is actually only a representation up to 
a sign factor. 

In two dimensions, the 5*0(2) group can be parametrized by a single 
parameter cf) 6 [0, 27r). For a given value of this parameter the rotation 
matrix is given by : 

«(#)=^ = ^) +e -n(#)=(_^ $») (27) 

where e a 6 is the standard antisymmetric matrix with e x 2 = 1. 

The corresponding matrix TZ(<f>) in the spinorial representation is TZ((f>) = 
ea " 12 *. In particular, if we set 7 1 = 03 and 7 2 = oi, where cij are the Pauli 
matrices, then a 12 = 02 and rotation matrix is : 

where e = io"2 is an antisymmetric tensor that is numerically identical with 
the one in fl27D . The difference, of course, is that the tensor in equation (|27|) 
has frame indices e a b whereas the one in (^HJ) has spinorial indices e£ . 

In order to fix the global sign of 7Z{4>), on should control the angle 4> in 
the range [0, 47r) rather than the usual [0, 2tt). This would require changing 
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Figure 2: Rotation of a local frame by 2tt. Even though the resulting frame 
configuration is obviously the same as before, spinor components can change 
their sign due to the sign ambiguity. 



continuously the angle and calculating the overall change / d<p keeping track 
of the number of 'full circles'. However, this cannot be done here since the 
relative angles between the frames Ci a are determined in the fundamental 
range [0, 2ti) only. 

The sign ambiguity also has topological consequences. Consider once 
more the regular, toroidal, flat lattice and choose on it a constant field of 
identical frames (see fig. Q). We first set Uij = t for all links. Trivially, if 
at a vertex % the frame is rotated by 2ir, the frame configuration does not 
change. However, because TZi(2ir) = — 1 in the spinorial representation, all 
links emerging from % acquire a negative sign Uji = — 1 according to the trans- 
formation law ( p6f) . The resulting 'sign field' is different from the original one 
but at the same time equivalent to it. By repeating this procedure in other 
vertices one can produce many different, but equivalent, sign configurations 
for the same field of frames. 

It is easy to see that a local rotation of a frame by 2ir preserves the 
overall sign of all elementary plaquettes, i. e. the product of signs of all links 
on the plaquette's perimeter. Thus, for any configuration obtained from the 
original one, all elementary plaquettes have a positive overall sign. We shall 
require this to be true in general, i. e. for any configuration of local frames 
on the lattice the sign of all elementary plaquettes is set to +1; this ensures 
that spinors remain unchanged by parallel transport around any elementary 
plaquette. This requirement is dictated by the underlying continuum theory, 
in which parallel transport of a spinor around a closed loop in a locally flat 
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Figure 3: A small deformation of a loop L (bold line) by an elementary loop 
P (dashed line), resulting in the loop V . 



patch leaves the spinor intact. Later on, for curved lattices, we shall modify 
this constraint so as to adjust it to the case where there is a deficit angle 
inside an elementary plaquette. 

Assuming that all elementary plaquettes have a positive sign we can prove 
now some simple topological theorems concerning the signs of loops on the 
lattice. 

It is convenient to define an auxiliary operation for loops on a lattice, to 
be called a small deformation of a loop. To deform a loop L, we pick an 
elementary plaquette P which shares at least one common link with L, and 
substitute the intersection L H P by the complementary part of P, resulting 
in a new loop L' = LUP-LnP (see fig. |).Q 

As with elementary plaquettes, we can define the overall sign of a loop 
as the product of signs of all links on the loop. One easily checks that the 
sign of the deformed loop U is the same as that of L - namely, the addition 
of P to L cannot change the sign because P has a positive sign by default, 
and the removal of the intersection L n P cannot change the sign because 
each link is 'removed twice' (once from P and once from L), so that the total 
number of removed links is always even. 

Any contractible loop can be obtained from the elementary loop by a 
sequence of small deformations. Thus all contractible loops have positive 
signs. 

1 Somewhat more precisely, we also have to require that the intersection L n P be 
connected, so as to avoid situations in which a small deformation splits a loop into two or 
more parts. 
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Figure 4: A non-contractible loop on a toroidal lattice with a constant frame. 
The single links drawn as bold lines all have transition matrices Uji = — 1, 
whereas all other links have Uji = 1; as a consequence, the loop has a negative 
overall sign. 



This is not, however, the case with non-contractible loops, which can take 
either sign. An example of a loop with negative sign is shown in fig. |] : if 
we choose Uji = — 1 for one complete row of links on the lattice (as in the 
figure) and Uji = 1 everywhere else, then any loop that encircles the lattice 
in the y direction passes through exactly one link with negative sign, and 
thus has a negative overall sign 0. 

Obviously, two sign configurations are equivalent if one can transform one 
into the other by a sequence of local rotations lZi(2ii) — — 1. Because local 
rotations do not change the sign of any loop, a configuration with at least 
one loop of negative sign cannot be equivalent to a configuration that has 
only loops of positive sign. In other words, the two sign configurations are 
topologically distinct. 

Now, using small deformations we can easily prove that all non- 
contractible loops encircling the torus in the same direction must have the 
same sign. This means, for example, that it is sufficient to calculate the sign 
of just one 'vertical' loop (which encircles the lattice in the y direction) to 
know the sign of all other vertical loops. More generally, the sign of a loop is 
not a property of a single loop but rather of all loops in the same homotopy 
class, i. e. those that can be obtained from each other by a sequence of small 

2 More generally, if a loop which encircles the lattice in the y direction goes back and 
forth having a sort of S shape, it may cross links with negative signs more than once. The 
number of crossings is however odd. 
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Figure 5: (Left) A lattice with toroidal boundary conditions. (Right) A lat- 
tice with the boundary conditions of a Klein bottle. The arrows indicate the 
directions in which the opposite edges are to be taken when joined together. 



deformations. On the torus there are two independent non-trivial homotopy 
classes of loops ('vertical' and 'horizontal') and, therefore, four distinct pos- 
sible sign configurations. These, in turn, correspond to four distinct spin 
structures. 

The statement can be generalized by observing that there are 2g indepen- 
dent classes of non-contractible loops on a surface with genus g, which means 
that there are 2 2g different sign configurations and thus the same number of 
spin structures. In particular, a lattice with spherical topology admits only 
one spin structure. 

On the other hand, on a non-orientable lattice one cannot globally define 
a field of orientable frames. An example of such a lattice is the so-called 
one-sided torus or Klein bottle, which is constructed in the same way as the 
standard torus but has different boundary conditions, as shown in fig. ||. It 
is possible to show that a frame transported along a closed path would have 
changed its handedness after a complete tour around the lattice. Because 
there does not exists a field of orientable frames, one cannot in this case 
define a spin structure or a Dirac operator. 

4 Local frames on a random lattice 

The form flTU| ) of the Wilson action is particularly simple not only because of 
the simple topology of the torus, which allows for the definition of a global 
frame, but also because of the regular geometry of the lattice which every- 
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where repeats the same simple motif. On an irregular lattice, local angles and 
link lengths change from point to point. This must be reflected in the con- 
struction of the hopping term, which depends on these local details through 
the co variant derivative. 

To make the geometrical part of the discussion as simple as possible, 
and to minimize the number of local degrees of freedom of the lattice, we 
restrict the discussion to equilateral random triangulations. This greatly 
reduces the number of local degrees of freedom, making the discussion more 
transparent and allowing us to focus on the interesting topological part of 
the problem. Let us mention, however, that the presented construction can 
be easily generalized to the case of variable link lengths and angles. 

On an equilateral triangulation, the local geometry is completely encoded 
in the connectivity of the lattice; all other details are fixed by the simple 
geometry of the equilateral triangle. In particular, the deficit angle at a 
vertex i is determined solely by its order : Aj = (6 — q^ir/ft. 

The local curvature of the lattice is concentrated in the vertices of the 
triangulation. The geometry becomes singular in these points and therefore 
it is difficult to provide a unique definition of a tangent space at the vertices. 
It is more convenient to define tangent spaces at the dual points of the lattice, 
i. e. at the centers of the triangles. Inside each triangle the geometry is locally 
flat and thus naturally spans a tangent space. We therefore locate all local 
frames, and also all fermionic fields, at the centers of the triangles. Each 
point i where a field is defined has then three neighbors, each of which at the 
same distance from i. The vectors pointing to the neighbors are also equally 
spaced in the angular variable, i. e. they are separated by angles 2%/3. 

Before defining the fermionic fields, however, let us discuss the properties 
of the field of oriented orthonormal local frames on such a random triangu- 
lation. An example of a triangulation decorated with frames is shown in fig. 

At each triangle % live two orthonormal vectors and such that 
^ia'^ib = Sab. Apart from the internal product there is also an external one A, 
which enables one to choose frames with the same handedness e^iAe^ > for 
all triangles. Now consider two neighboring triangles % and j, each endowed 
with its own frame and ej. The interiors of the two triangles together form 
a flat patch of the triangulation. One can think of the two frames as being 
two alternative frames for the same patch. One can calculate components 
of our objects in either one of them, and easily recalculate them when going 
from one to the other. To this purpose introduce 5*0(2) transition matrices 
Uij and Uji such that : 
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Figure 6: A small piece of a random triangulation with local frames. Ujk is 
the transition matrix between the frames at k and j, and qi is the order of 
the vertex i. 

One can repeat the same calculation for any pair of neighboring triangles 
and use it to transport a frame between any two points i\ and i n along an 
open path C = (zi, i 2 , ■ ■ ■ , i n ) : 

e i„ = Ui n i n _ 1 . . . Ui 3 i 2 Ui<ii x £% x = U(C) ti x . (30) 

Since we study a theory whose content is independent of the choice of frames, 
we are interested in the pertinent transformation laws and in quantities in- 
variant under local 50(2) rotations of the frames : e« — > e'j = i^ej. The 
object U(Cji) = Ujk ■ ■ ■ U n i for any open path between i and j transforms 
as : 

U{Cfl) - U'(C Jt ) = R.UiC^Rr 1 , (31) 

as one can see from (p9|) . In particular, for a closed path Lj beginning and 
ending at the same triangle i, U(Li) transforms as 

U(Lt) -> U'(Li) = RtUiL^Rr 1 , (32) 

and hence Tr U(Li) is an invariant. Moreover, this invariant does not depend 
on the choice of the initial point i of the loop, and is thus a property of the 
loop L itself. It is a geometrical quantity related simply to the total angle 
/ da by which a tangent vector is rotated when transported along the loop. 
On a flat lattice, this angle is a multiple of 2n. On a curved lattice the 
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situation is somewhat more complicated. In particular, for an elementary 
loop L q surrounding a vertex of order q, the loop invariant is 

-Tr U (L q ) = cos — = cos — — ^— = cos A q (33) 
2 3 3 

and contains information about the deficit angle A q , or equivalently about 
the curvature at the vertex. There are various possibilities to prove this 
statement; the proof outlined here offers us an opportunity to introduce an 
auxiliary construction which will be useful throughout the remaining part of 
the paper, especially when we shall lift the spin connection to the spinorial 
representation. 

Recall that the information about the local geometry of the lattice is 
stored in the form of three local unit vectors nfl pointing from i to its three 
nearest neighbors. There is, however, another and for the problem at hand 
more suitable way of achieving the same goal. Instead of the vectors 

themselves, one can equivalently consider the rotations that connect nf- to 
Ci. To introduce the rotation matrices, we first associate an entire frame 
with each of the three nearest neighbor vectors, treating n?} as the first basis 
vector of each corresponding frame. The second base vector of the frame 
is then automatically determined by the orthonormality condition. Now we 
have three particular frames njf jQ = (n^ l5 n^ 2 ) f° r the three neighbors j of i. 
The frames tv?} can be obtained from the local frame by a rotation Bj : 

nf = Bfei . (34) 

We refer to them as to the basic rotations at i. 

Now, it is convenient to decompose the connection matrices Uji into basic 
rotations Bj at i and B± at j. Letting them act first on the frame e^, one 

obtains the frame n?} . One then flip it to the frame using a rotation 
by 7r, which is represented by the matrix F = e e7T . Finally, using the inverse 
basic rotation at j one rotates it to ej. In other words, the transition from ej 
to ej (and vice versa) can be done in the following three steps (see fig. 0) : 

ej = \B?Y l FBf et, e, = [Bf)- X FB? e 3 . (35) 

Comparison with (|29| ) leads to : 

U,j = \BfY l FB? , U 3l = [B^FBf . (36) 

One can use this decomposition to calculate the loop invariants Tr U (L) : 

n 

Tr U(L) = TrY[U lk+llk =Trl[T lk (37) 

k=l k 
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Figure 7: A patch of two neighboring triangles, and the three nearest neigh- 
bor vectors for each of them. The same information can be provided by a 
rotation matrix between and the first basis vector e^i, shown as the flag 
emerging from the center of each triangle. In this example, the basic rota- 
tion B\ of frame j to the direction of its neighbor % is a rotation by 5ir/3, 
whereas the basic rotation of frame i to the direction of its neighbor j 
is a rotation by tt. 

where TT is an ordered product that runs through all vertices on the loop 
L = (hih, ■ ■ ■ ,i n ) with the cyclic boundary condition z n+1 = i\ and the 
rotation matrices 

T h = Bf k) \Bf k) = e (±) ^f e (38) 

correspond to the turn taken by the path at the triangle ik [|19|]. It depends 
on the turn-angle, which can be either +tt/3 if the path turns to the left 
or — 7r/3 if it turns to the right. In fact, on a equilateral triangulation, the 
sign (±)i fc determines completely the turn matrix T ik at the triangle ik- It 
does not depend on the particular orientation of the frame, because under 
rotation of the frame ik the basic rotations transform as : 

B (i k ) B^R ih , [B^]- 1 -> R^[B {ik) Y l (39) 

thus leaving the combination [_B(* fc )] _1 in T ik intact. 

An elementary loop around a vertex of order q turns exactly q times in 
the same direction. Thus we have 

^Tr U(L q ) = Tr e±*« = cos f = cos . (40) 

as claimed in (|3"B|). 
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5 The spinorial representation 



The next step is to lift the connections to the spinorial representation, 
Uij — > Uij. We continue to use the convention of denoting all rotation ma- 
trices in the spinorial representation by calligraphic letters : U — > U for 
connections, B — ► B for basic frame rotations, T — > T for turns and F — > T 
for flips. 

The starting point of the construction is the decomposition (|36]). If we 
write it in the spinorial representation, each matrix that occurs in this equa- 
tion is determined only up to a sign : 

e e<f> _^ ± e ^/2 (Hj)_ The idea is now 
to affix the spinorial representation of all matrices on the right hand side of 



(|36 ) with a positive sign : 



B = e^ -> £ = e e * /2 (41) 
F = e e7r = 1 -> ^ = e e7r/2 = e , (42) 

and keep the sign Sjj = ±1 as a separate variable for each link : 

U. . U .. = s .. [ B f]~h B? , Uji - ^ = [B? 5 ]-^ Bf . (43) 

We demand that parallel transport of a spinor along a given link and back 
does not change the spinor. We see that this is indeed the case, i. e. we have 
UjiUij = 1 if 

Sji s tj = -1 . (44) 

Using a similar calculation as the one which led to ( |40D one finds that in the 
spinorial representation the loop invariant for an elementary loop around a 
vertex is 

±TrU(L q ) = S Lq -cos^. (45) 

where A q is the deficit angle, and 5*^^ is a sign ±. The factor one-half in 
the argument of the cosine follows from The total sign of the loop, 

denoted by S , ^ g , depends on the choice of signs s^- in (|43"D and has to be 
calculated. We require that the signs Sy are chosen in such a way that for 
each elementary loop the sign is positive : 

Sl 9 = 1 • (46) 

Note that for q = 6 this requirement is natural, because the plaquette is 
flat, A 6 = 0, and as discussed before for a flat patch the parallel transport 
should be trivial : U(Lq) = t. Thus indeed we should have Sl 6 = 1. Also for 
other g's the requirement can be motivated. The geometry of an elementary 
plaquette corresponds to the geometry of a flat cone, which has a singularity 
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a 




Figure 8: The internal geometry of a set of triangles around a vertex is 
the same as that around the peak of a cone : it is flat everywhere except 
for a single point where the curvature is concentrated in a singularity. We 
can determine the sign of any loop around the cone if we first regularize this 
singularity by 'flattening' the cone, and find S = +1. 



at the peak. The elementary loop encircles this singularity at some distance 
r from the peak. One can regularize the singularity by smoothing the peak, 
i. e. replacing it by a differentiable surface (see fig. ||). 

In doing so, one deforms only a very small region within a distance of 
e around the peak, where e r. Now imagine that we shrink the loop, 
continuously decreasing its radius. Then TrW(r) and A(r) both change con- 
tinuously with r. In the limit r — > 0, the loop ends up on the top of the 
regularized part of the geometry which is flat. Thus, again S = +1 in the 
limit of r — > 0. This already is sufficient to have positive sign for all values 
of r, because in the course of continuous changing, the deficit angle A was 
changing continuously and hence the sign S could not have jumped between 
negative to positive values without making U discontinuous. In other words, 
S must keep the value +1 for all r. 

Because the regularized zone can be made arbitrarily small, we assume 
that the triangulated lattice, which corresponds to the limit e — > 0, inherits 
the property of the regularized geometry: the sign of any elementary loop is 
S Lq = +1 for any q. 

In order to enforce the constraint Si q = +1 for each plaquette, one has 
to establish a relation between Si q and the signs of links Sji. In analogy to 
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(pr?D , one can calculate the loop invariant in the spinorial representation as : 

n 

TrU(L) = Tr l[U ik+lik = ' Tr IK ■ ( 47 ) 

fc=i fc & 

Comparing this to the result pertinent for the fundamental representation 
(|37|), one finds that an additional product of link signs appears, as expected. 
But there is also another source of signs hidden in (f47|). It has its origin in 
the spinorial representation of the turn matrices T — > T. Surprisingly, and 
in contrast to the fundamental representation, the product of basic rotations 
depends on the position of the frame. More precisely, calculating the rotation 
corresponding to the turn taken by the path at if. one gets an additional sign 

% k = B^IB^T = z lk e^>^ (48) 

which was not present in the fundamental representation. 

The reason for the appearance of these new signs is the following: In the 
spinorial representation, the basic rotations are given by 

n. . — ph^'k+i^k 6 n. . — ph&k-iik 6 (aq) 

where 4> ik+lik and 4> ik _ lik are the angles between (e ifcl , n ik+li J and 
(ei k i,n ik _ lik ), respectively. Therefore, we have 

j: _ e i(& fc+ iife- < ^-i i fe +7r ) e = e K A ^ fc + 7r ) e . (50) 

By construction, <fii k+1 i k and </>i fc _ 1 i fc both lie in the range [0, 2tt). However, 
the difference A(j) ik = 4>i k+1 i k — 4>i k -iik can ne outside this range. In general, 
one has A0 ifc + 7r = ±n/3 modulo 2tt, but 2tc can be disregarded since 
e 2ire _ j n spinorial representation, however, due to the factor 1/2 one 
has (A0 ifc + 7r)/2 = ±7r/6 modulo ir, and this 7r cannot be ignored because 
e we = ±1. 

One has to calculate the exponents in (|50"D exactly and to find all possible 
values of A0, fe . There are six different cases, collected in fig. ^. 

The flag in each drawing represents the position of the vector e^i, with 
respect to which the angles are calculated. We call it the z-flag. For example, 
in the drawing (a) one has 4> ik+lik G [0,2^/3) and <f> ik _ lik = <pi k+1 i k + 4vr/3, 
which yields A0 ifc = — 47r/3 and thus the rotation matrix : 

7J k = e 3(-^/3+^) = e -f e . (5i) 

In the drawing (b) one has 4> ik+lik G [2tt/3, 47r/4) and <f> ik _ lik = <pi k+1 i k -2vr/3, 
so that A0j fc = 27r/3 and the rotation matrix is 

%k = e ^/3 + ,)e = e f e = _ e -| e (52) 
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Figure 9: The six different possibilities for a path to cross a triangle with a 
marked z-flag, constructed from the two possible directions of the path (left 
turn or right turn) and the three possible directions of the flag. The sign of 
A0j fe is determined by whether or not the auxiliary line to the right of the 
path crosses the flag. 









(a) 


-4tt/3 


+e -7r/6e 


(b) 


+2tt/3 


_ e -vr/6e 


(c) 


+2tt/3 


_g— 7r/6e 


(d) 


-2tt/3 


_|_ e +7r/6e 


(e) 


+4tt/3 


_g+i"/6 e 


(f) 


-2tt/3 


_|_ e +i"/6e 



Table 1: The difference of angles A0 ifc and the turning matrix % k in the 
spinorial representation for the six cases shown in fig. |9|. 
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The results for all six cases (a - f) are given in table [I]. Inserting them into 
the formula for the loop invariant (|47|) one obtains : 

TrW(L) = II<W* Tr I[% = II ' Tr \{e^^ (53) 

k k k k 

where z ik is the sign of % k . Setting : 

S l = - l[s ik+lik z ik , (54) 
k 

one finds : 

TrW(L) = -S L - Tr He (±)l ^ e . (55) 

k 

The relation flo"4] ) between the loop sign Sl, the link signs s, and the z-signs 
can be represented graphically in a very intuitive way. The signs Zi k tell on 
which side of the path lives the z-flag. If one draws an auxiliary line, as in 
fig. |], that runs along the right-hand side of the main path, then the sign zi k 
can be determined geometrically by choosing Zi k = —1 if the auxiliary line 
crosses the z-flag and z ik = +1 otherwise. Similarly, one can introduce a field 
of flags associated with the oriented links, and choose the sign Sji = — 1(+1) 
when the respective s-flag is (is not) crossed when one is going from i to j. 
Because for any given link the auxiliary path crosses the s-flag when going in 
one direction but not in the other, this choice leads to s^Sy = —1 as required 
by fl44]). The total sign Sl of the loop L is now given by the number of flags 
that are crossed by the auxiliary path : 

S L = (-1) 1+Fl . (56) 

As on the regular lattice, one can use the concept of small deformations of 
loops to prove some topological theorems for the signs of the loops. The fact 
that each elementary loop has S = +1 implies that two loops L, V that can 
be transformed into each other by a small deformation always have the same 
sign, Sl = Sl', because a small deformation changes the number of flags 
crossed by the loop by an even number (see fig. |T0|) . 

Thus, we see that if all elementary loops on the lattice have positive signs, 
all contractible loops have positive signs Sl = +1, too. Likewise, one can 
show that all loops belonging to the same homotopy class have the same sign. 
In other words, all the topological theorems we found for the regular lattice 
hold for the triangulated one as well. The remaining thing is to check that on 
a given lattice an assignment of the link signs s^, ensuring the positivity of 
all elementary loops signs, does always exist. That it is so for any discretized 
orientable 2D manifold in [IU| , |Tl"| . 
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Figure 10: A small deformation of a loop on a triangulated lattice. 



6 The Dirac- Wilson operator 

We now have all what is needed to construct the fermionic action (|19|). We 
start by casting the formula (El) for the hopping operator 



Hlj ~ 2 



n)i ■ 7 



Uii (57) 



into a form that depends on the field of orthogonal frames through the basic 
rotations. One can use equation ( |43|) to decompose the matrix Uij : 

U ij = 8 ij [Bf]- 1 eB^ . (58) 

Likewise, we write the vector n£ in terms of the basic rotations. By defini- 
tion, the basic rotations at point i relate the direction en of the frame to the 
directions of the links between i and its neighbors j : 

e* = [Bf]- 1 ^ . (59) 

In the spinorial representation fll7|) one can write : 

41 ■ 7 = \B?rVBf , (60) 
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where 7 1 = en ■ 7 is the gamma matrix associated with the first direction of 
the frame. As mentioned before, the gamma matrices have the same numer- 
ical values 7 1 = 03, 7 2 = <j\ in each frame on the triangulation. Inserting 
everything into (|5TD we eventually obtain : 

W y = ay [Bf]-^[l-7 1 ]eB? ) , (61) 

which defines the hopping term in the Dirac-Wilson operator on the trian- 
gulated lattice. 

To calculate the basic rotations, one has to find on each triangle the three 
angles between e»i and the nearest neighbor vectors n^; denote them by <\k . 
Each is defined in the fundamental range of the rotation group, [0, 2tt). Since 
physical quantities cannot depend on the choice of the field of frames, we are 
free to make the most convenient choice. Hence, we assume that in each 
triangle the vector en points to one of the vertices. This implies that the 
angles <f>^ can take only one of the three possible values - tt/3, tx or 5n/3 - 
which in turn makes the basic rotation matrices very simple : 

m *f ( c (i) s {i) \ 
Bf = e^= % ■% , (62) 



where 



c)> = cos-^- = — ,0,-— , s)> =sin-^- = (63) 



for (f)j = 7r/3, 7r, 57r/3, respectively. Inserting this explicit form of the basic 
rotations into (|61| ) leads to an extremely simple formula because (1 — 7 x )/2 is 
a projection matrix, which with our choice of 7 1 has only one non-vanishing 
element. Hence : 



= s v %%) %)%) ■ ( 64 ) 



In this form the Dirac-Wilson operator is easy to implement. For each pair of 
neighboring triangles j and i we first find the sign Sji and the angles between 
the z-flag and the dual link ji and calculate the appropriate trigonometric 
functions. For example, assuming Sy = 1 for the link ji in fig. |7| we have 
<j>i = 57r/3, (f)j = 7T. Hence : 
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The Dirac- Wilson operator is built from blocks like the above one, for each 
pair of indices representing neighboring triangles, and from 2x2 unit matrices 
for each pair of identical indices. Defining the adjacency matrix for triangles 
as : 

. | 1 if i and j are neighbors 
3 | otherwise 

one can write the Dirac- Wilson operator as : 

V i3 = -KA ij H ij + 5 ij l. (67) 

What are the properties of the Dirac- Wilson operator in this form? Consider 
the charge conjugation transformation : 

V - ^ c = C^ T , V - = -VC~ l , (68) 
where the matrix C is unitary and fulfills the requirements : 

C- 1 7 T C = -7, C T = -C. (69) 
One can check that the hopping operator (^l|) transforms as : 

cn'jC • n,,. (70) 

In two dimensions we can choose the standard antisymmetric matrix e as 
the charge conjugation matrix, C — e. It is convenient to use two different 
versions of e, one with lower indices e a p and one with upper indices e a/3 , but 
with the same numerical values : 

s 12 = e 12 = l, e aj e^ = -6?. (71) 

One can treat e as a simplectic form to raise or lower the spinorial indices : 

(Vc)« = e^, m a = ^^ a - (72) 

We recall that in the explicit index notation, the components of the spinor ip 
are denoted by ip a and those of ip by ip a . Furthermore, in this notation one 
can write : 

V:f = e^[V lJ f 1 . (73) 

In the implicit index notation one has to distinguish between different cases, 
namely T> for mixed indices, eD for only upper indices, and T>e for only lower 
indices, by displaying explicitly the action of e. 

The fact that the hopping operator is constructed from a projector implies 
in particular, that : 

HijHji = , HijUjiHij = Hij . (74) 
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The consequence of the transformation law ( |70|) is that : 

eHije = -U% (75) 

and, furthermore, that : 

(e V ..f = - eVji . (76 ) 
In index-explicit notation, this last equation reads : 

V$ = -V?? , (77) 

which means that the matrix Dff is antisymmetric in the double indices 
I = (ia) and J = (j(3) : V„ = -V, H . 



7 Second-quantized theory 

Quantum field theory of free Dirac fermions in a curved geometrical back- 
ground represented by a triangulation T is defined by the partition function : 

Z T {K) = J IJdVid 2 ^e-^^«]^ = \ V \ . (78) 

i 

The propagator is : 

(VwC) = / II dV*d 2 ^ VvC e-^Ae = [V -^ . ( 79) 

It transforms under a local change of frames e, — > e! i = RiCi as follows : 

(W<> - (VCVC) = MM^VwVd)- (80) 

Let us further explore the consequences of the symmetry with respect to the 
charge conjugation that is encoded in the transformation law ([70]). Introduce 
two families of Major ana fermions : 

1 1 - 

<t>i = ~(i/>c + il>) , 01 = 2^ + ^), 

02 = 7^c-^), 02 = ^c-VO. (81) 

They are charge self-conjugate : 0i c = (f)% and 02c — 02- This means that the 
components of 0i are not independent, likewise for 2 - The components are 
related : 

0° = 0/3^- (82) 



27 



as can be seen from ([721). We skipped the family index 1, 2 in the last formula. 

It is convenient to express the Dirac- Wilson action in terms of the Majo- 
rana families 0i and 02- Indeed, using equation ( [TDD one finds that the two 
families decouple : 

Sty, = \ E *M>i ~ K H = Wi) + S{4> 2 ) (83) 

* (ij) 

where 

S(0) = \ E " ^ E Wo-^ • (84) 

(ij) 

The two actions Sty, if)) and 5(0) appear identical to each other, but they 
differ in the number of degrees of freedom; in the latter case, is uniquely 
determined by 0. By changing the variables in the integration measure of 
(|78|) one can rewrite the partition function as a product of two identical 
factors : 

Z T (K) = J nd 2 0Hd 2 2l e- 5 ^- 5 ^ = [Z T (K)f (85) 

i 

where Zr{K) is the partition function for a single Majorana family : 

i 

= J n d Vie~* lo2J ^ = Pfaff[eD] . (86) 

i 

Here, eD is the antisymmetric matrix (|77|), which implies that the square of 
the Pfaffian is equal to the determinant of eT>, which is in turn equal to the 
determinant of T>. We can calculate the partition function for the Majorana 
fermions using the hopping parameter expansion. This leads to a geometrical 
interpretation of the model, as will be seen in the next section. 



8 Fermionic loops 

To find the hopping parameter expansion of Zt{K) let us first split the 
integrand into two parts : 

Zr(K) = J J] (d 2 0* e-fr*) I] (l + Ko/H.jOj) (87) 

i (ij) 

The first part is a product of independent one-point integrations with an 
exponential measure, whereas the second is a product over all oriented links 
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that connect neighboring points. Since we know from equations (|70| ) and 
(75) that for Majorana fermions : 



bjHji4>i = (piHij(f>j , 



it is convenient to rewrite the product in ( |87j) as a product over non-oriented 
links (ij) : 

Zt{K) = J II ( d ^ e ~ ¥ ^) II + IKhHijh) . (89) 

To do this, we have to require that terms like (piTiijtpj^jTijifpi do not occur 
in the expansion. Actually, they vanish because of fl74|) . 

The only non- vanishing integrals relevant to our problem are : 

e"^.l = l (90) 

and / 

' d 2 e~3**. (0-0) = 1. (91) 



These rules are used to calculate the integral of each term in the expansion : 
11(1 • 2 Ko,H, ,<:>,) = 1 • -IKY^oiH.jOj 

{ij) {ij) 

+{2Kf fcHiih ■ 4>kH hl <t>i + ■ ■ ■ (92) 

mm 

Consider the quadratic term on the right hand side. If j = k then, according 
to (pID , the integration over <pj yields : 

(fii'Hijrf'j • <t>jHji(j)i = <f ) i('Kij'Kji) ( l ) i ■ (93) 

Otherwise, if j ^ k, the integral vanishes. In general, one observes that the 
contribution of a term in the expansion (|92|) is non- vanishing only when all 
neighboring fields <pj ■ <ftk belong to the same point. Integration of these terms 
over all fields gives : 

where all ji in the chain are different. For the final integration to yield 
something non- vanishing one must have j\ = j n . Finally : 

C(L) = -Tr Hj lj2 Hj 2 j a ■ ■ ■ Hj 7l _ ljl . (95) 
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This contribution can be graphically represented by a closed loop L = 
0i> 32, ■ ■ ■ ,jn-i,ji) of length n. On the other hand, integration over a field 
<pk associated with a vertex k that does not lie on any loop contributes a 
factor of 1 (13). 



In summary, all terms of the expansion that survive the integration (|89|) 
can be represented graphically as diagrams consisting of closed loops. These 
loops do not back-track or touch each other. A configuration consisting of / 
loops Li, L2, ■ ■ ■ , Li with total length n = rix + . . . + ni contributes a term 

(2K) n C(L 1 )C{L 2 )...C(L l ) (96) 

to the partition function. 

One can calculate the contribution C (L) of a single loop L in a way similar 
to that used to obtain the loop invariant (f47"D , i.e. by extracting the total 
sign of the loop and expressing the remaining product in terms of turns 
at the vertices (EH). The result is : 



C(L) = -Tr UH ih+llih = S L -Tr U% k hl - V) (97) 



k k z 



The difference between this expression and the one for the loop invariant (|53| ) 
is that now in addition to the turn matrix a projection operator appears in 
the product. Inserting the explicit form of the turn matrix Tj = e ±e?r//6 and 
of the projector (1 — 7 1 )/2 = (1 — o"3)/2, one obtains : 




C(L)=S L [^-\ . (9i 



This is again similar to the result found for the loop invariant (|55|), but 
with two differences. First, one now has instead of — Sl- Second, in 
the calculation of the loop invariant the turn angles enter the result with a 
sign ± depending on whether the path turns left or right, whereas here the 
projector leaves only the cosines of the rotation matrix, which depend on the 
absolute value of the turn angle. Thus, each turn contributes a factor +\^3/2 
independently of its direction. Since a loop makes a turn at each vertex, the 
number of turns in a loop is simply equal to the loop length, which gives 

Inserting this result into (p6[), one finds that the contribution of a loop 
configuration of total length n is : 

S total ■ (V3K) n , (99) 

where 

Stotal = U S L- (100) 
L 
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Figure 11: A configuration of fermionic loops. 



On a lattice with spherical topology all loops L have a positive sign Sl = 1 
and therefore S to tai = 1 for each loop configuration. 

On a torus, the sign of the contribution depends on the spin structure. 
Assuming periodic boundary conditions in both directions (++), all loops 
from any non-trivial homology classes, contractible or not, have Sl = 1, 
and again S to tai — 1 fo r any loop configuration. The standard notation 
is used here: the spin structure is referred to by the signs of independent 
classes of non-contractible loops. On the torus there are two classes and 
therefore four possibilities (ss f ), with s, s' = ±. Plus/minus corresponds 
to periodic/antiperiodic boundary condition for spinors transported along 
loops in this class. With anti-periodic boundary conditions in any direction 

- (H — ), ( — h), or ( ) - any non-contractible loop circling the lattice in 

this direction has a negative sign Sl = —1. Thus, all of these three cases 
can produce unwanted configurations with a negative contribution to the 
partition function. More generally, any configuration that has an odd number 
of non-contractible loops circling the lattice in an anti-periodic direction has 
a negative total sign Stotai = _ 1- 

Yet another possible choice of boundary conditions imposes summation 

over all spin structures - (++), (H — ), ( — h), and ( ) - in the partition 

function. This operation is called GSO projection, and in many cases seems 
to be the most physical choice. Negative contributions are not a problem 
in this configuration with an odd numbers of loops in one of the 

non-trivial homotopy classes, say in the first class of non-contractible loops, 

has St ota i = 1 for (++) and (H — ), but S total = -1 for ( — h) and ( ). The 

summation over all cases yields zero. More generally all 'bad' contributions 
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Figure 12: Domain walls versus loops on a torus. A non-contract ible loop 
on a torus, like for instance the upper curve in the figure, cannot be a part of 
the domain wall configuration of Ising spins unless there is a partner curve 
in the same class of loops in this configuration, like for example the lower 
one. In general, domain-wall configurations of 2D Ising model have an even 
number of loops in each non-trivial class of non-contractible loops. 



to the partition function cancel out in the GSO projection. 

Configurations with an odd number of non-contractible loops in at least 
one direction cannot correspond to Ising model domain wall configurations, 
because only an even number of these domain walls is crossed when one is 
performing a round trip on the lattice (see fig. 12). From this point of view, 
the loop cancellation in GSO projection is very physical. Before discussing 
this point in more detail, a more careful look at the properties of the loop 
signs is needed. 



9 The GSO projection 

As discussed in the preceding sections, the global properties of the Dirac- 
Wilson operator on a two-dimensional compact manifold are closely related 
to the signs of the fermionic loops. Self-consistency requires a positive sign 
for all elementary fermionic loops, and this in turn implies a positive sign 
for all contractible loops. Non-contractible loops, on the other hand, are 
not subject to this restriction. In fact, it is the ensemble of signs of all 
independent non-contractible loops that defines the spin structure of the 
manifold. 

In this section, it will be shown that the sign of any loop on the lattice is 
uniquely determined by the signs of a minimal number of independent non- 
contractible loops. Stated differently: the signs of all loops on the manifold 
are completely encoded in the manifold's spin structure. 
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So far, we discussed the loops without self-crossings only, for the simple 
reason that on a triangulation no other loops occur in the hopping parameter 
expansion of the Majorana-Dirac- Wilson fermions. On the other hand, we 
already encountered self-crossing implicitly in the calculation of the invariants 
TrW(C) (f|7|), since they can be defined on loops of any kind, including the 
self-crossing ones.Q 

For this reason, and also for the sake of completeness, we shall now discuss 
the signs and topological properties of loops in a general context, and restrict 
them to self-avoiding loops only when necessary. We require the sign of a 
loop to be a property of its homotopy, which means that we have to modify 
the definition of the sign fl56|) to : 



where Cll is the number of self-crossings of the loop L. Of course, for any 
contractible loop this must still result in a positive sign, independently of the 
number of self-crossings. A few examples of contractible loops with various 
numbers of self-crossings are shown in fig. [H| It is easy to verify that the 
auxiliary line running along the right hand side of the loop crosses an odd 
number of flags in the first two cases and an even number of flags in the last 
two. 

Let us now return to the operation that we called a 'small deformation'. 
So far, we have considered only deformations that do not induce self-crossings 
(see for example fig. [H]). These deformations will be called even. It is 
convenient to introduce also an odd version of a small deformation, where 
an elementary plaquette is again used to deform the loop, but with a a self- 
crossing like in fig. [14]. The two kinds of small deformations differ by the 
orientation of the plaquette that is used to deform the loop. Neither kind 
changes the overall sign of the deformed loop. 

One can introduce equivalence classes of loops that can be obtained from 
each other by a sequence of small deformations. A class of loops equivalent 
to a loop A will be denoted by [A]. Inside this class, [A] even denotes the 
sub-class of loops that can be obtained from [A] by a sequence of an even 
number of small deformations. 

Let us define the loop merging operation, that acts on a set of equivalence 
classes of loops. Take two loops A 6 [A] and B 6 [B], and deform both 
of them smoothly until they have a common link. If this common link has 

3 It is convenient to think of a self-crossing on a lattice not as a meeting at exactly one 
vertex, but rather as a sort of smeared overlapping that may occupy one or more links of 
the lattice. In particular, on a lattice with only vertices of order three, there are no exact 
one-vertex self-crossings; the most localized ones still occupy at least one link. 




(101) 
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Figure 13: Examples of contractible loops with and without self-crossings on 
a square lattice. The numbers of crossed flags and self-crossings are F\ — 9, 
C\ — for the first example; F 2 = 9, C 2 = for the second (the flag in the 
center of the figure is crossed twice!) ; F 3 = 8, C 3 = 1 for the third; and 
F3 = 12, C3 = 1 for the last one. The result is a positive sign in all cases. 
Note that in the last three examples, the flag at the vertex in the center, 
where four links of the loop meet, is crossed an even number of times by the 
auxiliary line. 




Figure 14: Even and odd versions of a 'small deformation'. In the upper 
figure, the two loops contain a common link of opposite orientation, causing 
the two versions of the link to 'cancel out' in the resulting deformed loop. 
In the lower figure, the common link has the same orientation in both loops, 
causing it to appear twice in the deformed loop and thus introducing a self- 
crossing. 
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Figure 15: The loop merging operation. Two loops A e [A] and B e [B] 
are smoothly deformed until they share a link. They are then joined either 
by erasing the link or by a self-crossing, depending on the link's relative 
orientation in both loops. The resulting loop belongs to a new class of loops 
[A- B). 



an opposite orientation in both loops, erase it and form a loop out of the 
remaining links. Otherwise, leave the link as it is and join A and B by a 
self-crossing (see fig. [T5p. The resulting loop belongs, by definition, to a new 
equivalence class of loops [A ■ B] . 

The product of loop classes defined in this way has a unity element in the 
class of contractible loops [E], for which [A ■ E] = [E ■ A] = [A]. Counting 
the number of crossed flags before and after the loop merging (steps 2 and 3 



in fig. |15|), one finds : 



F [a .b] m = 2 F [A] + F [B] + 1 + C AB , (102) 

where Cab is the number of crossings of the loops A and B. The equation 
implies the law of sign composition : 

S[A-B] = S[ A ] S[B] • (103) 



Indeed, in the upper drawing in fig. 15 the loop merging does not intro- 



duce any additional self-crossing, Cab = 0, and the number of crossed flags 
changes by 1 modulo 2, whereas in the lower figure one additional self-crossing 
appears, Cab = 1, and the number of crossed flags changes by modulo 2, 
i.e. it remaines unaltered. The factors coming from the flag count and from 
the number of additional self-crossings compensate each other, and the above 
simple composition law ( [103|) follows. 

Thus, the set of equivalence classes of loops forms a group with respect to 
the loop merging operation. On a two-dimensional compact manifold, this 
group contains a minimal set of independent classes of non-contractible loops 
[Hi], i = 1, . . . , 2g, where g is the genus of the manifold (see fig. |H)|). 
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Figure 16: Independent classes of non-contractible loops on a 2d manifold 
with genus g. 




Figure 17: Loop merging on a torus. Take two loops from the classes [Hi] 
and [H2] and smoothly deform them until they share a link, then merge them. 
The resulting loop [Hi ■ H2] circles the torus in both directions [Hi] and [H2] 
simultaneously. Its sign is the product Si ■ S2. 



This minimal set has the nice feature that all other classes can be created 
from [E], [Hi], and their inverses [-f^] -1 by use of the loop merging operation. 
In other words, one can decompose any loop in terms of [E] and [ifj], and 
then use equation (|103j) to calculate the sign of this loop as a product of 
signs of the Hi. 

Let us illustrate this with a few examples. For simplicity, denote the signs 
of the classes in the minimal set with Si = S^] = Sma-i. 

Consider first a loop which goes around a torus in two distinct homotopy 
directions simultaneously. Such a loop can be obtained by loop merging of 
the classes [Hi] and [H 2 ], as shown in fig. [T7[ Note that the loops shown 
in the figure do not self-cross; nor does the resulting loop Hi ■ H 2 . This 
might seem surprising at first, given that the loop merging itself introduces a 
crossing, Ch 1 h 2 — 1- But indeed one can see that the original loops, even if 
not self-crossing, do cross each other. In general, any loop from [Hi] always 
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Figure 18: Loop merging on a double torus. The sign of the loop in the lower 
figure can be calculated by observing that it can be created by a merging 
of the loops drawn in the upper figure. This can be done by first deforming 
the two loops until they have a common link, then erasing this link, and 
smoothly deforming the remaining loop. 



crosses any loop from [H 2 ] an odd number of times. In the resulting merged 
loop, these crossings become self-crossings, so that the product has an even 
number of self-crossings overall. This in turn means it can be deformed by 
a sequence of an even number of small deformations to a non-self-crossing 
loop. 

As a general definition, one can state that a class [A] crosses a class [B] 
if the number of crossings between any two representatives A and B is odd. 
By this definition, the classes [Hi] and [H 2 ] cross each other, as do any two 
of the classes [i^i-i] and [H 2 i] shown in fig. [16|. This concept will be useful 
in a while in the context of the GSO projection. 



Another example of loop merging is shown in fig. [TBI. The loop C in the 
lower drawing is obtained by merging Hi and if 3, so its sign can be calculated 
as the product Sc = S^S^- 

Let us apply the sign composition law to the calculation of the partition 
function of Majorana fermions (|89|) . In the hopping expansion, one generates 
non-self-crossing loops only. Denote the number of loops from a given class 
S[c] on a configuration by N\c] • Then the total sign of this configuration can 
be written as : 

Stotal = II (S[C\) ■ ( 104 ) 

[G] 
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After GSO projection : 

2 2 a 



Z gso = J_ j2 £((±)i.(±) 2 >->(±) 2g ) ; (105) 



{(±)l,...,2 ff } 

the total contribution of the configuration is proportional to 

w Gso ~n^( i +(- i ) Ar[c] ) • ( ioe ) 

[C] 2 

To see this, note first of all that the sum over the signs Si of the classes 
[Hi\ can be replaced by a sum over the signs S[c] of the classes [C] present 
in the configuration, since all these loops do not cross and are independent 
from each other. Summing over all signs S\c] means that each loop of each 
non-trivial class [C] occurs an equal number of times with plus and minus 
signs, which eventually leads to the last formula. In a sense, the action of 
the GSO projection factorizes into a product of independent actions for the 
loops of each non-trivial class on the configuration. 

The last equation also tells us that all configurations with an odd num- 
ber of loops from any non-trivial class have a vanishing contribution to the 
GSO projection. Physically, this means that the projection removes all loop 
configurations which cannot represent domain wall configurations. 



10 Topology of the Ising model 

We shall consider now the Ising model with nearest neighbor interactions, 
focusing on the issue of the exactness of the duality transformation between 
the model defined on a triangulation and on its dual graph, respectively, and 
emphasizing the topological aspect of the duality. Furthermore, the relation 
between the Ising and the fermionic model will be discussed. 

To distinguish between a triangulation and its dual, we attach a star to 
symbols referring to the triangulation, while the unstared symbols refer to 
the dual lattice. 

With this convention, the partition function of the Ising spins living on 
the triangulation reads : 

nr.Ga.) = 4t +) (&) = E e* , (107) 

to*} 

where <7j, = ±1 are spin variables located at the vertices z* of the triangula- 
tion. As we shall see later discussing boundary conditions for the Ising model, 
the partition function ( |107| ) corresponds to the partition function with the 
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Figure 19: Ising spins on the triangulated lattice, and the corresponding 
domain walls drawn as loops on the dual graph. 



spin structure (++). Therefore we additionally denoted it by in 
the last equation. 

Any spin configuration on the triangulation can be graphically repre- 
sented as a configuration of loops on the dual graph. Namely, for any link 
connecting two spin variables of opposite sign, a if = — , one can draw that 
link's dual as a part of a loop. It is easy to see that the result will be loops 
surrounding domains of aligned spins (fig. |19|). 

One can calculate the statistical weight of every loop configuration. For 
this purpose, it is convenient to rewrite the partition function as : 

n T .(ft) = e ^ 2Nf3 * J2 e^cW^-- 1 ) , (108) 
to*} 



which can be obtained from ( II 07] ) by subtracting unity from the link in- 
teraction energy. Since we consider triangulations without boundaries, the 
numbers of links and dual links are equal, Nl = Nl*, arid related to the 
number, N, of triangles by Nl = 3/2N. The subtraction of unity in each 
interaction term is compensated by adding an appropriate constant factor in 
front of the sum in (|108| ). The contribution to the sum of a term (a^a^ — 1) 
is if <7j„ = (jj, , and 2 if ^ cr^ . Thus, the sum in the exponent gives twice 



the number of domain wall links (denoted as bold links in fig. [19]), which is 



equal to the total length n of all loops on the configuration. Therefore : 

n T ^*) = 2e i / 2Nf3 *Y. e ~ w * n i ( 109 ) 

{L} 
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where the sum runs over all loop configurations on the dual graph (which are 
identical to the loop configurations of the fermionic model discussed in the 
previous section). The additional factor of 2 in front of the sum reflects the 
fact that each loop configuration represents two distinct spin configurations 
which can be obtained from each other by a simultaneous flip of all spins 

For a non-spherical topology, some attention has to be paid to non- 
contractible loops. Consider once more a toroidal triangulation. A con- 
figuration with an odd number of non-contractible loops does not form a 
domain wall configuration of the Ising model and therefore does not appear 
in (|109|) . The same is true of the fermionic model if we perform the GSO 
projection. Therefore the equivalence between the models is exact : 

Z^°(K) = 2e~V 2N ^-n T M), (HO) 

if we set 

JlK = e" 2/3 * , (111) 

as can be seen by comparing ( |99|) and ( |109| ). This statement holds for an arbi- 
trary triangulation of a two-dimensional orientable manifold without bound- 
ary. 

Consider now the Ising model with spins <7j living on the vertices of the 
dual lattice, or equivalently at the centers of the triangles of the original 
manifold (in other words, the spins are located at the same spots as the 
Majorana fields fa discussed before). The partition function reads now : 

T (/3) = ^/Ecy)^'. (H2) 
{<*} 

Performing the strong coupling expansion leads to the formula : 

Ortf) = cosh(/3) 3 / 27V/3 E II (1 + tanh(/5)) , (113) 

in analogy to the hopping parameter expansion ( j39f) in the Majorana field 
theory. The integration rules for Ising spins : 

\ y L 1 = 1 > ^5> = °> l$> 2 = i ; (114) 

CT=± cr=± cr==b 

are completely analogous to those for the fermions fl90D, fl9T|)-H 

4 One would see a difference with the fermion rules on a lattice with vertex orders greater 
than three, because then one could also have terms like 1/2 J2a=± 0-4 = whereas the 
corresponding terms in the Majorana model are zero, / d 2 (j) e"5<^ • (fxfxfxj) = 0. However, 
in our case the order of the dual lattice vertices is three by construction. 
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Thus, calculating the strong coupling expansion, one again finds a sum 
over the same loop configurations : 

n T ([3) = (2cosh(/3)) 3/2JV/3 ^T(tanh/3) n . (115) 

{L}' 

More precisely, for a spherical lattice the loop configurations occurring in 
this sum are identical to the domain wall configurations of the Ising model 
defined on the triangulation. However, this is not true for topologies of higher 
genus, where configurations with an odd number of non-contractible loops 
from the same homotopy class occur in the strong coupling expansion ( |115| ). 
This is why we have put a prime on the sum, to distinguish the set of these 
configurations from the set of domain walls ( |109| ). If we again take the torus 
as an example, we see that the sum in (|115|) also contains configurations 
with a single loop, or with an even number of loops circling the torus in the 
indirection. This kind of loop configuration is also produced in the hopping 
expansion of the fermionic model if we restrict it to the spin structure with 
periodic boundary conditions. Thus, in this case we have : 

Z { + +) (K) = (2cosh((3)y 3/2NI3 ■ n T ((3) (116) 

if we set 

V3K = tanh(/3) . (117) 

The equivalence also holds for topologies of higher genus if we choose this 
spin structure for the fermionic model. 

As expected, both Ising models are almost dual to each other. The only 
difference comes from topological contributions related to non-contractible 
loops. In fact, one can make the two models exactly equivalent by a sort 
of a 'GSO projection' for the Ising field. Contrary to the projection in the 
fermionic model, which appears as a natural option because the model has 
several possible spin structures, its introduction here is somewhat artificial. 

Again, take the torus as an example. Originally, we have only one version 
of the Ising model, which corresponds to the spin structure (++) ( |107| ). Now, 
we attempt to define a model that can reproduce the three other structures. 
Let us start with the spin structure ( — h), corresponding to an anti-periodic 
boundary condition in the first homotopy direction. On this lattice, choose a 
non-contractible loop circling the torus once in the second homotopy direction 
(see fig. |20| ). We call this an anti-ferromagnetic line. All links (ij) that 



cross this line will be called anti-ferromagnetic and denoted by AH 
other links will be called ferromagnetic and denoted by + - We define the 
partition function as follows : 



ft$r%) = V A E «>+ CT -V . (118) 
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Figure 20: For the dual Ising model on the torus, define an anti-ferromagnetic 
line in the indirection as a non-contractible loop circling the torus in 
this direction. The Ising interaction for a given link is defined as anti- 
ferromagnetic or ferromagnetic depending on whether or not it crosses the 
anti-ferromagnetic line. 



In the strong coupling expansion, each ferromagnetic link contributes a factor 
+ tanh(/3), and each anti-ferromagnetic link, a factor — tanh(/3). Each non- 
contractible loop in the if i-direction has an odd number of anti-ferromagnetic 
links, so its contribution will be — tanh n (/3), whereas each contractible loop 
and each non-contractible loop in the ^-direction has an even number 
of anti-ferromagnetic links, thus contributing + tanh™(/5). In other words, 
this prescription gives exactly the same sign factors as those occurring for 
fermionic loops on the torus with spin structure ( — h). In the same manner, 
one can also introduce an anti-ferromagnetic line in the ifi-direction, to pro- 
duce a model corresponding to a (H — ) spin structure. Finally, a model with 

an anti-ferromagnetic line in both directions gives us a ( ) spin structure. 

Summing over all four cases, one obtains a model with a partition func- 
tion : 



n 



GSO 
T 



(P) = - + + 4~ +) (£) + 4" _) (/5)) , (H9) 



which is exactly dual to the Ising model Qqt + \(3*) that has its spin variables 



defined on the vertices of the triangulation ( |107| ), and is equivalent to the 
model of Majorana fermions with CSO-projection. 

For a lattice size going towards infinity, the difference between Q^ + ((3) 
and Qj> SO (/3) becomes negligible. As already explained, the difference comes 
only from the non-contractible loops. These loops can be regarded as hav- 
ing a one- dimensional entropy, in the sense that they can be ordered by a 
one-dimensional index that represents their position on the lattice. Because 
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X 

Figure 21: Fermions on the honeycomb lattice. 



of this, they become less and less important when the system size grows. 
Therefore, in the thermodynamic limit one expects an exact duality between 
( |107| ) and ( |112j ) even without extending the model to the spin structures 

(H — ), ( — h), and ( ). We introduced this extension here to ensure exact 

duality, i. e. a one-to-one map, between the two models even for systems of 
finite size. 

Generalization of this Ising model 'GSO construction' to higher genus 
topologies is straightforward. In order to simulate a spin structure with 



an antiperiodic boundary in a given direction Hi (see fig. W), one simply 
introduces an anti- ferromagnetic line in the direction Hj that crosses Hi. 
Altogether, this creates 2 2g different spin structures. 



11 Two examples 

As a first example, consider the Dirac- Wilson action on a regular triangula- 



tion of the two-dimensional plane (fig. |21]). The fermions live at the triangle 



centers, on a regular hexagonal lattice. Because the lattice is flat, we can 
choose a global frame, i. e. with the same directions e\ and e 2 at each vertex. 
To fix the signs, we also choose the flag assignments, which can likewise be 
done in a translationally invariant way. 

One can easily write down the fermionic action for this model. Choose 
an elementary cell as in fig.[H]. It consists of two distinct sites : A and B. 
The lattice can be constructed by shifting the elementary cell by multiples 
iid\ + 22^2 of the fundamental shift vectors di = n + ni, d 2 = n + n 2 
constructed from the the link vectors : 

n = (0,1), ni=(^,|), and n 2 = ±) . (120) 
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The components of the vectors n\ and 712 are expressed in the global frame 
(X, Y) shown in the figure. The position of the cell is referred to by the 
double integer index i = (z 1; i 2 ). With this notation the action is written as : 

K 2 

S = -J2 [^i+dA 1 + n d ■ l)i>i,B + AM 1 - n d ■ l)lpi+d,A 

1 i d=l 

K r - 



+ ^Y1 [VMVM + ?Pi,Blpi,B 



:i2D 



Since the plane is non-compact, topological effects are not relevant. From the 
discussion in the previous sections we know that for Majorana fermions the 
model with this action is equivalent to the Ising model with spin variables 
living at the vertices and at temperature f3* given by (|1 1 1|) , and likewise to the 



Ising model with spins at the centers of the triangles and at temperature (3 
given by ( |1 1 7| ) . The critical temperature corresponds to the critical hopping 



parameter, for which the fermions become massless. This critical value is 
easily found to be : 

K CT = ^ , (122) 

because each vertex on the dual lattice, where the fermions are living, has 
three neighbors. Thus, the critical temperatures for the Ising models is : 

& C r = -^hi^5, / 9 cr = iln(V3 + 2), (123) 



in agreement with the known results ]23| . 

A second example we want to discuss shortly here is the discretization of 
the Majorana field coupled to two-dimensional gravity. It is well-known that 
the integration measure over the metric field on a two-dimensional manifold 
can be represented by a sum over all equilateral triangulations. If we dress 
each triangulation in this sum with the fermion field, we effectively obtain a 
theory of Majorana fermions coupled to two-dimensional gravity. This theory 
is given by the partition function : 

Z{K) = Y J Z T {K), (124) 

T 

with the sum running over all triangulations with a fixed topology. For non- 
spherical lattices, one should sum in addition over spin structures. 

We can use now the equivalence between the Majorana-Dirac- Wilson ac- 
tion and the Ising model to substitute, triangulation by triangulation, all 
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Figure 22: Comparison of the results for the energy density of Ising field 
computed from MC simulations of the Ising model (line) and of the cor- 
responding qunatity eq. ( |128| ) from MC simulations of the fermionic model 
(crosses). The error bars are smaller than the symbols used. 



terms Zt in the sum. We again obtain an exact map between the Ising 
model and the model of fermions coupled to gravity. The Ising model, how- 



ever, is exactly solvable |]24" |; in particular, the critical temperature is [25 



Pcr = l^, & cr = iln^, (125) 

which means that the Majorana fermions are massless when the hopping 
parameter K is : 

K cr = ^e~^ = ^. (126) 
393 K J 

This, again, is an exact result. The equivalence of the two models opens the 
possibility of studying numerically the properties of the Dirac- Wilson oper- 
ator coupled to gravity. In fact, one can use the Ising model as a generator 
for triangulations, and then dress the configurations with local frames and 
z and s flags to calculate T> on each of them. Since the Dirac- Wilson oper- 
ator depends on the triangulation, one gets a model of dynamical fermions 
interacting with the fluctuating geometry. Using the Ising model as a Monte 
Carlo generator for configurations is many orders of magnitude more efficient 
than a generator referring directly to the fermionic action, since using the 
latter requires calculating the Pfaffian flSB]) in each single Monte Carlo step, 
an extremely costly operation in terms of CPU time. One can easily con- 
vince oneself by simulating small systems that the two generators do indeed 
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produce the same results but differ enormously in algorithm efficiency. In 
fig. ^ we compare the average energy of the Ising field calculated in the two 
different ways : (a) directly using the Ising model : 

<J hi^T, (127) 




Nd(3* 

i 

or (b) using the equivalence ( |11(J| , |111| ) : 



3 1 dK dlnZ 1 / 1 _ A 

e *~~2~NW^K~ " ~2~\2N\ K ) t (128) 

where A a are eigenvalues of the Dirac-Wilson operator T> on the given tri- 
angulation T. In the derivation of the last formula we made use of the 
relations : 

dK \ dK / T 2\ dK j T K ) 

The two methods yield the same results. Using the trick with the Ising 
model as a generator of triangulations one can extend the MC simulations 
to larger systems in order to investigate the properties of the spectrum of 
the Dirac-Wilson operator on dynamical triangulations. The results of these 
investigations has been presented elsewhere [p6fl . Here let us only quote 
a result for the finite size scaling of the pseudocritical hopping parameter 
K* defined as the value of the hopping parameter for which a mass gap is 
minimal. By the mass gap we mean the center of mass of the distribution of 
the smallest positive eigenvalue of the Majorana-Dirac- Wilson operator eD. 
The numerical results can be well fitted to the finite size scaling formula : 

K, = K 00 + ^ (130) 

where K^ = 0.3756(16), and k = 1.03(30), a = —0.9(5). The parameter 
K^ corresponds to the critical value of the hopping parameter in the ther- 
modynamic limit. As one can see it agrees with the theoretical prediction 
K cr = 0.3746... given by the equation ( |126| ). 



12 Conclusion 

The topological properties of a fermion field on discretized two-dimensional 
compact manifolds were discussed at length. The exact equivalence between 
the model of Majorana- Wilson fermions and the Ising model was established. 
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An exact duality relation for the Ising model on a compact manifold was also 
found. 

It would be important to generalize the construction to higher-dimensio- 
nal simplicial manifolds. Having done this, one would then be able to attack 
the problem of quantum gravity interacting with a fermionic field. So far, 
it has only been possible to couple integer spin fields to four-dimensional 
simplicial gravity jl6[ p|. Such a theory is known to have problems with 
the continuum limit [ 17 1 , which could reflect the fact that higher-dimensional 
gravity does not exist without a proper cocktail of matter fields coupled to it. 
If this were true, the addition of fermions might perhaps help solving these 
problems. 

It is straightforward to generalize parts of the construction presented in 
this paper to higher dimensions. In particular, one can associate with each 
four-dimensional simplex an orthonormal oriented frame and basic rotations, 
and out of them one can easily build the transition matrices and spin connec- 
tions. However, the problem of lifting this construction to the half-integer 
representation leads to additional complications. 

One of the reasons is that the topological problem is by itself more com- 
plicated in four dimensions. The question of whether a manifold admits a 
spin structure, which is equivalent to the question of whether it is possible to 
define globally a Dirac operator on it, is in general related to the existence 



of a non-trivial second Stiffel- Witney form |2l|. For two- and three- 
dimensional manifolds, this reduces to the orientability question. In four 
dimensions, however, there are manifolds, like for example the projective 
space CP(C 2 ), which are orientable but possess a non-trivial Stiffel- Witney 
form, and thus do not admit any spin structure. In an attempt of extending 
our construction to a higher dimensional manifold not admitting any spin 
structure, the topological obstruction would manifest itself as the impossi- 
bility to adjust the local degrees of freedom so as to assign positive signs to 
all elementary plaquettes. 
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